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Abstract 

Scaled momentum distributions for the strange hadrons Kg and A/A were mea- 
sured in deep inelastic ep scattering with the ZEUS detector at HERA using an 
integrated luminosity of 330 pb -1 . The evolution of these distributions with the 
photon virtuality, Q 2 , was studied in the kinematic region 10 < Q 2 < 40000 GeV 2 
and 0.001 < x < 0.75, where x is the Bjorken scaling variable. Clear scaling vio- 
lations are observed. Predictions based on different approaches to fragmentation 
were compared to the measurements. Leading-logarithm parton-shower Monte 
Carlo calculations interfaced to the Lund string fragmentation model describe 
the data reasonably well in the whole range measured. Next-to-leading-order 
QCD calculations based on fragmentation functions, FFs, extracted from e + e~ 
data alone, fail to describe the measurements. The calculations based on FFs 
extracted from a global analysis including e + e~, ep and pp data give an improved 
description. The measurements presented in this paper have the potential to fur- 
ther constrain the FFs of quarks, anti-quarks and gluons yielding K s and A/ A 
strange hadrons. 
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1 Introduction 



The jet fragmentation and hadronisation processes through which coloured partons become 
bound in colour-neutral hadrons cannot be described within the framework of perturbative 
QCD (pQCD). Several approaches have been developed which attempt to build a bridge between 
the fixed-order partonic cross sections and the observed hadrons. Two of the most successful 
and widely used approaches are the Lund string model [1] and the fragmentation functions 
(FFs) [2-6]. The Lund string model, relying on a large number of parameters, is interfaced to 
leading-logarithm parton-shower Monte Carlo models. The FFs are parameterisations of the 
hadronisation process within the standard framework of leading-twist collinear QCD factorisa- 
tion, in a similar way to that of the parton distribution functions (PDFs), and are convoluted 
with the predicted partonic cross sections. 

Extensive studies of the fragmentation properties of the hadronic final state have been performed 
in e + e~ [7-19], pp [20-23], pp [24] and deep inelastic ep scattering 1 (DIS) [25-33] data and have 
provided information about the fragmentation and hadronisation processes. The measurements 
provided tests of pQCD and showed that scaling violations are observed. In addition, the 
comparison of the measurements in different reactions indicated an approximately universal 
behaviour of quark fragmentation. 

In a previous publication [28], the ZEUS Collaboration presented high-precision measurements of 
inclusive charged-hadron production. Next-to-leading-order (NLO) QCD calculations, based on 
different FFs obtained from fits [34-36] to e + e~ data, from fits [37] to e + e~ , pp and pp data and 
from fits [38,39] to e + e~, pp and ep data, were compared to the measurements. The predictions 
based on the different FFs are similar and fail to provide a good description of the measurements 
over the full range of applicability of the calculations. The parameterisations [38,40,41] of the 
FFs for strange hadrons, such as Kg and A, are so far largely unconstrained. The ep data 
presented in this paper have the potential to constrain these FFs over a wide kinematic range. 

In this paper, the scaled momentum distributions for Kg and A hadrons 2 are presented for the 
first time in DIS. The scaled momentum is defined = 2P Breit / v / Q2, where P Breit is the 

particle momentum in the Breit frame and Q 2 is the photon virtuality. The Breit frame [42,43] 
is the frame in which the exchanged virtual boson is purely space-like, with 3-momentum q = 
(0, 0, — Q), providing a maximal separation between the products of the beam fragmentation 
and the hard interaction. The measurements were performed in the current region of the Breit 
frame, which is equivalent to one hemisphere in e + e~ annihilations, as functions of Q 2 and x p . 
Next-to- leading-order predictions, based on different FFs, and leading- logarithm parton-shower 
Monte Carlo calculations, interfaced with the Lund string fragmentation model, were compared 
to the measurements. 



1 Here and in the following, the term "electron" and the symbol "e" denote generically both the electron 
(e~) and the positron (e + ), unless otherwise stated. 

2 Here and in the following, the notation A includes both the particle and its antiparticle unless otherwise 
stated. 
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2 Theoretical framework 



In lowest-order QCD, three processes contribute to the DIS cross section, namely the Born 
(V*q —> q, with V* = 7*, Z*), the boson-gluon-fusion (V*g — > qq) and QCD-Compton-scattering 
(V*q — > qg) processes. The cross section for the production of an observed hadron, H, in the 
final state in DIS can be expressed in QCD, using the factorisation theorem, as 

a{ep ->• e + H + X) = ^ f j/p (x,Q)®a jj/ (x,Q,z)®F H/j ,(z,Q), 

j,j'=q,q,g 

where the sum runs over all possible initial (final)-state partons j (j 1 ), fj/ p are the proton PDFs, 
which give the probability of finding a parton j with momentum fraction x in the proton, &jj> is 
the partonic cross section, which includes the matrix elements for the three processes mentioned 
above, and Fjjui are the FFs, which give the probability that a hadron H with momentum 
fraction z originates from parton j'. The scaled momentum variable x p is an estimator of z. As 
for the PDFs, the FFs include contributions from quark, anti-quark and gluon fragmentation. 
Absolute predictions for the FFs cannot be calculated; however, the dependence of the FFs on 
the scale Q is calculable in pQCD and governed by renormalisation group equations, similar as 
for the PDFs. 

The range of applicability of the FFs is limited to medium to large values of z, since the 
assumption of massless hadrons leads to a strong singular behaviour for z — > 0. At small z, 
finite mass corrections are important. However, the inclusion of small-z mass corrections is not 
compatible with the factorisation theorem and thus the FFs with mass corrections cannot be used 
with fixed-order calculations. A possible solution is to introduce a posteriori mass-correction 
factors to take this effect into account [37]. 

A large improvement in the precision of the ingredients of the calculations has been achieved 
in the last few years. Matrix elements up to NLO accuracy are available for many processes; 
for DIS, this corresponds to 0{a 2 s ). Parton distribution functions have become increasingly 
more precise, largely due to the high-precision HERA data. On the other hand, FFs, though 
increasing in accuracy [34-41], still lack the precision of the proton PDFs. 

The data most widely used to extract the FFs comes from e + e~ annihilations into charged 
hadrons [7-19]. These data are very precise and the predicted cross sections do not depend 
on PDFs. However, they do not provide information on how to disentangle quark and anti- 
quark contributions to the FFs and the gluon fragmentation remains largely unconstrained. In 
addition, the e + e~ data have poor statistics at large z, leading to large uncertainties in this 
region of phase space. Several parameterisations of the FFs exist [34-36]. 

In the last few years, new one-particle inclusive measurements coming from both pp collisions [20- 
23] and DIS [44] became available. The inclusion of these data in the extraction of the FFs yields 
a much more complete picture of the fragmentation process and provides a direct handle on 
quark, anti-quark and gluon contributions. A global QCD analysis of e + e~, pp and DIS data is 
now available for several hadrons [38,39] . This global FF set agrees with the previous extractions, 
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based on e + e data alone, in the regions of phase space which are also well constrained by e + e 
data alone. 



3 Experimental set-up 

A detailed description of the ZEUS detector can be found elsewhere [45,46]. A brief outline of 
the components most relevant for this analysis is given below. 

Charged particles were tracked in the central tracking detector (CTD) [47-49], the microvertex 
detector (MVD) [50] and the straw tube tracker (STT) [51]. The CTD and MVD operated in 
a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD consisted 
of 72 cylindrical drift-chamber layers, organised in nine superlayers covering the polar-angle 3 
region 15° < 9 < 164°. 

The MVD silicon tracker consisted of a barrel (BMVD) and a forward (FMVD) section. The 
BMVD contained three layers and provided polar-angle coverage for tracks from 30° to 150°. 
The four-layer FMVD extended the polar-angle coverage in the forward region to 7°. After 
alignment, the single-hit resolution of the MVD was 24 /im. The transverse distance of closest 
approach (DCA) to the nominal vertex in XY was measured to have a resolution, averaged over 
the azimuthal angle, of (46 © 122/py) /im, with pt in GeV. The STT covered the polar-angle 
region 5° < 9 < 25°. For CTD-MVD tracks that pass through all nine CTD superlayers, the 
momentum resolution was ct(pt)/pt = 0.0029£>t © 0.0081 © 0.0012/p^, with pt in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [52-55] covered 99.7% of the to- 
tal solid angle and consisted of three parts: the forward (FCAL), the barrel (BCAL) and the 
rear (RCAL) calorimeters. Each part was subdivided transversely into towers and longitu- 
dinally into one electromagnetic section (EMC) and either one (in RCAL) or two (in BCAL 
and FCAL) hadronic sections (HAC). The smallest subdivision of the calorimeter was called 
a cell. Under test-beam conditions, the CAL single-particle relative energy resolutions were 
a{E)/E = 0.18/\/£ for electrons and <r{E)/E = 0.35/v^E for hadrons, with E in GeV. 

The energy of the scattered electron was corrected for energy loss in the material between the 
interaction point and the calorimeter using the small-angle rear tracking detector [56, 57] and 
the presampler [56,58]. 

The luminosity was measured using the Bethe-Heitler reaction ep — > ejp by the luminosity 
detector [59-61] which consisted of two independent systems. In the first system, the photons 
were detected by a lead-scintillator calorimeter placed in the HERA tunnel 107 m from the 
interaction point in the lepton-beam direction. The second system was a magnetic spectrom- 



3 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing towards the 
centre of HERA. The coordinate origin is at the nominal interaction point. 
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eter arrangement [62], which measured electron-positron pairs from converted photons. The 
fractional uncertainty on the measured luminosity was 1.8%. 

4 Event selection 

The data used in this analysis were collected during the running period 2005-2007, when HERA 
operated with protons of energy E p = 920 GeV and electrons of energy E e = 27.5 GeV, and 
correspond to an integrated luminosity of 330 pb" 1 . The criteria to select DIS events are 
described below. 

A three-level trigger system [46, 63] was used to select events online. It relied on the presence 
of an energy deposition in the CAL compatible with that of a scattered electron. At the third 
level, an identified electron [64] with an energy larger than 4 GeV was required. 

Offline, the kinematic variables Q 2 , inelasticity, y, and the Bjorken scaling variable, x, as well as 
the boost vector to the Breit frame were reconstructed using the double-angle (DA) method [65], 
which uses the angles of the scattered electron and of the hadronic system. 

Deep inelastic scattering events were selected by the following requirements: 

• E' e > 10 GeV, where E' e is the scattered-electron energy; this ensures a reconstruction 
efficiency above 95% and a purity of the scattered electron of ~ 100%; 

• Ue < 0.95, where y e is the inelasticity estimated from the energy and angle of the scattered 
electron; this excludes events with spurious electrons in the forward region, which are pro- 
duced predominantly by photoproduction; 

• yjB > 0.04, where yjB is the inelasticity estimated using the Jacquet-Blondel method [66]; 
this rejects events for which the DA method gives a poor reconstruction; 

• 35 < 5 < 60 GeV, where 5 = J2( E i ~ p z % ) and E { is the energy of the i-th CAL cell, P Zi 
is the momentum along the Z axis and the sum runs over all CAL cells; this removes the 
phase space where photoproduction background and events with initial-state radiation are 
expected; 

• \Z vtx \ < 50 cm, where Z vtx is the Z component of the position of the primary vertex; this 
reduces background from events not originating from ep collisions; 

• \X\ > 12 and \Y\ > 12 cm, where X and Y are the impact positions of the scattered electron 
on the RCAL, to avoid the low-acceptance region adjacent to the rear beampipe; 

• the analysis was restricted to events with 10 < Q 2 < 40000 GeV 2 and 0.001 < x < 0.75. 
These requirements selected a sample of 2.16 • 10 DIS data events. 



4 



5 Kc and A selection and reconstruction 



The strange hadrons K s and A were identified via the charged-decay channels, K s — > tt + it~ 
and A — > pir~ (A — > pir + ). The candidates were reconstructed using two oppositely charged 
tracks associated with a displaced secondary vertex. In the case of the Kg, the mass of the pion 
was assigned to both tracks. For the A, the mass of the proton was assigned to the track with 
the largest momentum, whereas the mass of the pion was assigned to the other track, since the 
proton always has a larger momentum than the pion for A baryons with momentum larger than 
0.3 GeV. 

All tracks were required to be in the region of high CTD acceptance, |ry track | < 1.75, where 
r] = — ln(tan#/2) is the pseudorapidity in the laboratory frame and 6 is the polar angle with 
respect to the proton beam direction. The tracks had to pass through at least three CTD 
superlayers and were required to have transverse momenta P^ rack > 150 MeV. 

The analysis was restricted to the current region of the Breit frame by boosting the tracks to 
this frame and requiring p^ rcit < 0, where pjp clt is the longitudinal momentum of the track 
in the Breit frame. The combined four- vector momentum of the two tracks in the Breit frame, 
pBreit^ wag use j t reconstruct x p . 

Additional selection criteria, similar to those used in a previous analysis [67], were applied to 
the selected candidates to maximise the purity of the sample with a minimum loss of statistics. 
These requirements were: 

• dca < 2 cm, where dca is the distance of closest approach of the two tracks forming the 
candidate; 

• x 2 /dof < 5 for the x 2 °f the secondary vertex fit; 

• M(e + e~) > 60 MeV, to eliminate background from photon conversion; 

• M(pvr) > 1121 MeV (M(7r+7r~) < 475 MeV), to eliminate A (K s ) background from the K s 
(A) sample; 

• &2D < 0.03 rad, where O2D is the collinearity angle in the XY plane between the K s (A)- 
candidate momentum vector and the vector defined by the interaction point and the K$ (A) 
decay vertex; 

• &3D < 0.04 rad, where O^d is the collinearity angle between the Kg (A)-candidate three- 
momentum vector and the vector defined by the interaction point and the K s (A) decay 
vertex; 

• Lxy > 0.5 (1) cm, where Lxy is the distance between the K s (A)-candidate decay vertex 
and the primary vertex in the transverse plane; 

• P^ A > (<) 0.11 GeV, where Pj? A is the projection of the pion momentum onto a plane 
perpendicular to the Kg (A) momentum direction (the Podolanski- Armenteros variable [68] ) . 
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Figures 1 and 2 show the dca, #2D> &3D an d Lxy distributions for data and Monte Carlo (see 
Section 6) for Kg and A candidates, respectively. The description of the data by the Monte 
Carlo simulation is adequate. 

Figure 3 shows the M(ir + Tr~) and M (jm) distributions after these requirements. A small amount 
of background is observed. The fit shown in Fig. 3 is for illustration only. The number of Kg 
(A) candidates in each bin of x p and Q 2 was estimated by counting the entries in the signal 
region, 472 — 522 (1107.0 — 1124.5) MeV, and subtracting the number of expected background 
entries. The latter was determined from a linear fit to the sideband regions 403 — 422 and 
572 - 597 (1086.0 - 1098.2 and 1133.2 - 1144.4) MeV, also indicated in Fig. 3. There were 
806 505 (165 875) Kg (A) candidates in the data sample. In the current region of the Breit 
frame, there were 238 153 K s and 40 728 A candidates. A Monte Carlo study showed that 6% 
of the selected A candidates come from higher-baryon decays. 

6 Monte Carlo simulation 

Samples of Monte Carlo (MC) events were produced to determine the response of the detector 
and to correct the data to the hadron level. The MC samples were also used to compute 
predictions to be compared to the measurements. 

The generated events were passed through the Geant 3.21-based [69] ZEUS detector- and 
trigger-simulation programs [46]. They were reconstructed and analysed by the same program 
chain as used for the data. Particles with lifetime longer than 3 • 10~ n s, such as K s and A, 
were treated as stable at generator level and their decays were simulated by Geant. 

Neutral current DIS events were generated using the program Lepto 6.5.1 [70]. Radiative effects 
were estimated using the Heracles 4.6.6 [71,72] program with the Djangoh 1.6 [73,74] inter- 
face to Lepto. Heracles includes QED corrections for initial- and final-state radiation, vertex 
and propagator terms, and two-boson exchange. The QCD cascade was simulated using the 
colour-dipole model (CDM) [75-78], including the leading-order QCD diagrams as implemented 
in Ariadne 4.12 [79,80] and, alternatively, with the MEPS model of Lepto. Fragmentation into 
hadrons was performed using the Lund string model [1], as implemented in Jetset 7.41 [81-84]. 
The default parameter setting from the DELPHI/EMC [19,85] tune was used for the hadroni- 
sation. The CTEQ5D [86] proton PDFs were used for these simulations. 

7 Corrections and systematic uncertainties 

The measured scaled momentum distributions were corrected to the hadron level and to the QED 
Born level. The correction factors were calculated bin-by-bin using the MC samples described in 
Section 6. The correction factors take into account: (i) the event-selection efficiency for the cuts 
listed in Section 4, but for the Q 2 and x requirements; (ii) the efficiency to identify the Kg and 
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A decays, as specified in Section 5; (Hi) the migrations between bins due to detector resolution, 
which affects in particular the transformation to the Breit frame; (iv) the relevant branching 
ratios; and (v) the extrapolation to the full phase space. The factors calculated in the measured 
(x p ,Q 2 ) bins varied from 0.05 (0.05) to 0.18 (0.11) for K% (A) candidates, and reached « 0.25 
for candidates with momentum in the range 1 — 1.5 GeV and — 1 < 77 < 1; the lowest values were 
found for high Q 2 and x p values. Bins with an acceptance smaller than 0.05 were not used in 
the analysis. The QED correction factors were computed using the Monte Carlo samples; they 
arc below 5% for Q 2 < 100 GeV 2 and increase to a maximum of 20% at the highest values of 
Q 2 . 

The total systematic uncertainties on the scaled momentum distributions are larger than the 
statistical uncertainties in most bins. The statistical uncertainties themselves vary significantly 
over the kinematic range. For K$ (A), they are at the 1 (4)% level at low Q 2 and between 10 
to 90% (20 to 70%) over the x p range at large Q 2 . Many of the systematic uncertainties were 
observed to scale with the statistical uncertainty. In the following list, typical values of the 
uncertainties on the scaled momentum distribution are given separately for Kg and A, either as 
percentages of the statistical uncertainty or as absolute values: 

• imperfections in the simulation causing uncertainties on DIS event reconstruction and se- 
lection resulted in uncertainties of ^39% and ^49% of the statistical uncertainties. This was 
evaluated by modifying the selection cuts within the experimental resolutions. At low Q 2 , 
the variation of the cut on yjB from 0.04 to 0.07 resulted in large uncertainties exceeding 
these typical values; 

• an uncertainty of —2% in the overall tracking efficiency resulted in absolute uncertainties of 
+4% and +4%; 

• detector-alignment uncertainties affecting the calculation of the boost vector to the Breit 
frame resulted in uncertainties of ^25% an d 1^5% of the statistical uncertainties. This was 
evaluated by varying separately the simulated polar angle of the scattered electron and of 
the hadrons by ±2 mrad; 

• uncertainties on the K$ and A selection efficiency resulted in uncertainties of 7Jgo% and lgo% 
of the statistical uncertainties. This was evaluated by varying the cuts listed in Section 5: 
the dominant effects were due to modifications of the cuts on Q2D to 0.015 and 0.06 and #3D 
to 0.02 and 0.08; 

• assumptions concerning the details of the simulation of the hadronic final state resulted 
in absolute uncertainties of 1 3 % and 1^5%- At large Q 2 , these uncertainties were larger 
and exceeded TJgo% and 1^5%. This was estimated by using MEPS instead of CDM in the 
calculation of the correction factors; 

• background-subtraction uncertainties resulted in absolute uncertainties of and 1|%. At 
large Q 2 , the uncertainties exceeded these typical values and were as high as ±35% for both 
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K° s and A. This was evaluated by varying the size of the background window by ±40% and 
changing the background fit function from first to second order. 



The systematic uncertainties were added in quadrature for each bin. The total systematic 
uncertainty is dominated by the uncertainty in the simulation of the hadronic final state. At 
low Q 2 , the overall tracking efficiency also contributes significantly. At high Q 2 , the uncertainties 
related to the Kg and A selection are important. 



8 NLO QCD calculations 

Next-to-leading-order QCD calculations, which combine the full NLO matrix elements with 
the proton PDFs and FFs as explained in Section 2, were compared to the measurements. 
For the comparison, the observable x p is assumed to be equal to the variable z. For each 
bin in x p and Q 2 , a prediction was derived by numerical integration over the multiplicities 
d 2 m(H)/dzdQ 2 , with m(H) the number of H per DIS event. Two sets of calculations based on 
different parameterisations of the FFs were used. The first set was obtained from fits to e + e~ 
data and based on the program Cyclops [87], called "AKK+Cyclops" [36,37]. The second 
set was obtained from a global fit to e + e~ , pp and ep data, called "DSS" [38]. It was used only 
for K$ predictions. 

The AKK+Cyclops calculations were performed using Q as the factorisation and renormal- 
isation scales; the number of active quark flavours was set to re/ = 5; the proton PDFs were 
parameterised using the CTEQ6M sets [88] and Aqcd was set to 226 MeV. The calculations 
were done assuming massless particles. Hadron-mass effects [89] for Kg and A were included 
as correction factors [37]. The influence on the shapes of the calculated scaled momentum dis- 
tributions due to the mass effects is expected at small values of x p and Q 2 , as explained in 
Section 2. 

In the DSS calculations, the scaled momentum distributions were obtained by convoluting the 
NLO DSS set of FFs together with the MRST NLO [90] PDFs and appropriate NLO coefficient 
functions. For these calculations, A'^-mass corrections were not included. The predictions were 
computed as ratios for each bin, such that a later combination of bins is not possible [91]. 

The uncertainty from terms beyond NLO was estimated by varying the renormalisation scale 
by factors 0.5 and 2. The uncertainties from FFs could not be evaluated so far; it is to a 
certain extent represented by the differences in the predictions of AKK+CYCLOPS and DSS. 
In addition, it should be noted that the DSS FFs were extracted from data at low Q 2 and that 
the fits are thus almost unconstrained at high Q 2 [38]. 
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9 Results 



Scaled momentum distributions, (l/N)(n(H)/ Ax p ), with n(H) the number of H (K s or A), N 
the number of DIS events in a given Q 2 bin and Ax p the width of the x p bin, were measured in 
the current region of the Breit frame. The distributions are presented as functions of Q 2 and x p 
in the kinematic region of 10 < Q 2 < 40000 GeV 2 and 0.001 < x < 0.75. 

Figure 4 shows the scaled momentum distributions for K s as functions of Q 2 in different regions 
of x p . The results are also presented in Table 1. The data show clear scaling violation. This 
behaviour is expected on the basis of the QCD description of the parton evolution with increasing 
Q: the phase space for soft gluon radiation increases, leading to a rise of the number of soft 
particles with small x p . 

The predictions from the CDM and MEPS models, based on leading-logarithmic matrix elements 
plus parton shower and the Lund fragmentation model, as described in Section 6, are compared 
to the measurements in Fig. 4. They describe the shapes of the distributions fairly well while 
overestimating the overall production of K s by 10 to 20%. 

The NLO QCD calculations, based on full NLO matrix elements and the fragmentation-function 
approach described in Sections 2 and 8, are also compared to the measurements in Fig. 4 for 
x p > 0.1. For z < 0.1, the calculations become singular. 

The AKK+Cyclops calculations, based on FFs extracted from e + e~ data alone, fail to describe 
the measurements. These calculations predict a much too high K s rate but for x p > 0.6. These 
discrepancies might come from the fact that the FFs used in these predictions have a poorly 
constrained gluon contribution, which is dominant at low x p . 

The DSS calculations, based on FFs extracted from a global analysis, give a good description 
of the measurements for x v > 0.3 and 10 < Q 2 < 40000 GeV 2 . The prediction for this region 
of phase space is mainly constrained by pp data, which sufficiently constrain the FFs at high 
x p . At lower x p , the DSS calculations fail to describe the data. This can be explained by the 
fact that the DSS fit in this region of phase space is mostly unconstrained by the available data. 
Thus, the measurements presented in this paper will help to improve significantly such global 
fits in this region of phase space. 

Figure 5 and Table 2 show the scaled momentum distributions for K s as functions of x p in two 
regions of Q 2 . The predictions of CDM and MEPS give a good description of the data. In 
both regions of Q 2 , both NLO calculations predict too-steep spectra. At low Q 2 , this effect is 
especially pronounced. 

Figures 6 and 7 show the scaled momentum distributions for A. The results are also presented in 
Tables 3 and 4. Scaling violations are clearly observed. The predictions of CDM and MEPS give 
a reasonable description of the measurements, but overestimate the overall A rate by ~ 20%. 
The AKK+Cyclops NLO calculations fail to describe the measurements. As seen in Fig. 7, 
the predicted spectra in x p are, as in the case of K s , significantly too steep. 
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ZEUS has previously published measurements of scaled momentum distributions for inclusive 
charged particles in DIS [28]. These measurements are dominated by the contribution from 
charged pions. Figure 8 shows the scaled momentum distributions presented in this paper 
together with those from the inclusive charged particles analysis in the kinematic region of 
0.1 < x p < 0.4 as functions of Q 2 . For Q 2 > 100 GeV 2 , all distributions show a plateau. At 
lower Q 2 , and especially at low x p , sizeable mass effects are expected. This is clearly visible. 
For 0.1 < x p < 0.2, the value of (l/N)(n(H) / Ax p ) drops to 10 (20)% of its maximum value 
for A {Kg), while for inclusive charged particles, the (1 / N)(n(H) / Ax p ) value is still 40% of the 
plateau value at the lowest Q 2 accessible. 



10 Summary and conclusions 

Scaled momentum distributions for K s and A hadrons were measured for the first time in ep DIS. 
The distributions were measured in the Q 2 range from 10 to 40000 GeV 2 and 0.001 < x < 0.75. 
Scaling violations were clearly observed for both the Kg and A hadrons. 

Next-to-leading-order QCD calculations, based on different parameterisations of the FFs, were 
compared to the measurements. The predictions based on FFs extracted from e + e~ data alone 
fail to describe the measurements. Those predictions based on a global analysis which include 
e + e~, pp and ep data give an improved description of the measurements. However, they predict 
a too high production rate of K s and A hadrons at low x p and Q 2 . The measurements presented 
in this paper have the potential to constrain significantly the FFs for the strange hadrons K s 
and A. 
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Q 2 (GeV 2 ) 


0.0 < x p < 0.1 


0.1 < x p < 0.2 


0.2 < x p < 0.3 


10 - 40 


0.031±0.00lt[];$] 2 


0125±0.00lj#{jo3 


0. 144±0.001±g;8?§ 


40 - 160 


0.171±0.002l{j;{£jj 


0.392±0.003iHJo 


0.283±0.002l°;°oe 


160 - 640 


0.551±0.010+°;° 2 ^ 


0.612±0.010±°;° 2 ? 


0.306±0.007l°;°o9 


640 - 2560 


1.141±0.038ig:ra? 


0.618±0.030^;^ 


0.309±0.029l^n 


2560 - 10240 


1.878±0.168lg;^ 


0.834d=0.217 r ±gl2re 


0.115±0.062±g;gH 


Q 2 (GeV 2 ) 


0.3 < x p < 0.4 


0.4 < x p < 0.6 


0.6 < x p < 1.0 


10 - 40 


0.1112±0.0008±°;°°^ 


0.0130±0.0004+°;°°^ 


0.0132±0.000ll°;°°°J 


40 - 160 


0.1571±0.0019^| 


0.0671 ±0.0009tg;g^ 


0.0109±0.0003lg;^ 


160 - 640 


0.1585±0.0060iHo5? 


0.0548±0.0027lg:8oi4 


0.0073±0.0008lHo?? 


640 - 2560 


0.1053±0.0217tg;g2i4 


0.0558±0.0141+°;°^ 


0.0029±0.0022ig;gg^ 



Table 1: The measured scaled momentum distributions (l/N)(n(Kg)/Ax p ) as 
functions of Q 2 in different regions of x p . The statistical and systematic uncertain- 
ties are also shown. 
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Xp 


10 < Q 2 < 100 GeV 2 


100 < Q 2 < 40000 GeV 2 


0.0 - 0.1 


0.0488±0.0006lg;8oi2 


0.4841±0.0063j£oii7 


0.1 - 0.2 


0.1618±0.0009lg;gg^ 


0.5740±0.006llH?32 


0.2 - 0.3 


0.1648±0.0009+°S 


0.3140±0.0048+°;°^ 2 


0.3 - 0.4 


0.1183±0.0007t{S 


0.1588±0.0037lg;gJ^ 3 


0.4- 0.5 


0.0751±0.0006±g;ggfi 


0.0760±0.0027lHoi7 


0.5 - 0.6 


0.0452±0.0004+°;°°^ 


0.0408±0.0022ig;S8i! 


0.6 - 0.7 


0.0260±0.0003+°;°°^ 


0.0182±0.0015+°;°°^ 


0.7- 0.8 


0.0150=b0.0002±g;ggg| 


0.0101=h0.0012±g;gg2| 


0.8 - 0.9 


0.0073±0.000llg:oooi 


0.0034±0.0007lg:ooo! 


0.9 - 1.0 


.0032±0 .000 1 ±S;8Sge 


0.0020±0.0006t{j;{j{j{j? 



Table 2: The measured scaled momentum distributions (l/N)(n(Kg)/Ax p ) as 
functions of x p in different regions of Q 2 . Other details as in the caption to Table 1. 
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Q 2 (GeV 2 ) 


0.0 < x p < 0.1 


0.1 < x p < 0.2 


0.2 < x p < 0.3 


10 - 40 


0.0025±0.0002j£ooo2 


0.0122±0.0004j£ooo8 


0.0189±0.0005lHou 


40 - 160 


0.0189±0.0014±°;°°^ 


0.0650±0.0021±g;gg|g 


0.0656±0.0018ig;g^ 


160 - 640 


0.0995±0.0095±g:88ra 


0.1666±0.0099±g;g^j 


0.0960±0.0066l°;°o38 


640 - 2560 


0.2313±0.0427lg;g567 


0.1966±0.0414l°;^ 


0.1038±0.0346lg;^ 


2560 - 10240 


0.7416±0.4386tg;|^ 


0.1962±0.1268t[J;?395 




Q 2 (GeV 2 ) 


0.3 < x p < 0.4 


0.4 < x p < 0.6 


0.6 < x p < 1.0 


10 - 40 


0.0184±0.0004+°;°^ 


0.0106±0.0002+°;°°^ 


0.0021±O.OOOllo™ 


40 - 160 


0.0453±0.0015lg;g^? 


0.0198±0.0007tg;^S 


0.0020±0.0002l°;°°°? 


160 - 640 


0.0641±0.0074lg;g^ 


0.0139±0.0025lg:ooli 


0.0041±0.00091Htos 


640 - 2560 


0.0653±0.0436lg;{j^ 







Table 3: The measured scaled momentum distributions (l/N)(n(A)/Ax p ) as 
functions ofQ 2 in different regions of x p . Other details as in the caption to Table. 1. 



17 



Xp 


10 < Q 2 < 100 GeV 2 


100 < Q 2 < 40000 GeV 2 


0.0 - 0.1 


0.00437±0.00034±8:oooi7 


0.08307±0.00638j#oorae 


0.1 - 0.2 


0.01822±0.00059i[j;H 


0.13904±0.00466lg;E 


0.2 - 0.3 


0.02456±0.00059+°™ 


0.09489±0.00275+°;™ 


0.3 - 0.4 


0.02173±0.00049lg;^ 


0.056013±0.00188lHo353 


0.4 - 0.5 


0.01387±0.00035±g;ggg§g 


0.02950±0.00135±g;gg|fS 


0.5 - 0.6 


0.00913±0.00026+°;™ 


0.014640±0.00095±°;™ 


0.6 - 0.7 


0.00483±0.00018+°;™ 


0.00534±0.00051+°;°°°34 


0.7-0.8 


0.00245±0.0001ltg;S 


0.00178±0.00028l^oo49 


0.8 - 0.9 


0.00096±0.00006±8;g8^ 9 


0.00056±0.00019lHoois 


0.9 - 1.0 


0.00038±0.00004+°;™ 


0.00032±0.00013^;S 



Table 4: The measured scaled momentum distributions (l/N)(n(A)/Ax p ) as 
functions of x p in different regions ofQ 2 . Other details as in the caption to Table. 1. 
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Figure 1: The normalised (a) dca, (b) # 2 d, (c) 9zd and (d) Lxy distributions 
for data ( dots ) and Monte Carlo (histograms ) for K s candidates. 
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Figure 2: The normalised (a) dca, (b) 62D, (c) &3D an d (d) Lxy distributions 
for data (dots) and Monte Carlo (histograms) for A candidates. 
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Figure 3: (a) The 7r + 7r~ invariant-mass distribution for K° s candidates (dots), 
(b) The pn" /p7i + invariant-mass distribution for A/ 'A candidates (dots). In both 
(a) and (b), the solid line represents an indicative fit by two Gaussians and a (a) 
linear and (b) quadratic background function. The solid vertical lines indicate the 
signal window used in the analysis. The dashed lines indicate the two sideband 
regions used for the background subtraction in each kinematic bin. 



21 



ZEUS 




Figure 4: The measured scaled momentum distributions (l/N)(n(Kg)/Ax p ) as 
functions of Q 2 in different regions of x p (dots). The inner error bars represent 
the statistical uncertainty. The outer error bars show the statistical and systematic 
uncertainties added in quadrature. In some bins, the error bars on the data points 
are smaller than the marker size and are therefore not visible. For comparison, the 
NLO predictions of AKK+Cyclops (dark-shaded band) and DSS (light-shaded 
band) are also presented. The bands represent the theoretical uncertainty. The 
predictions from CDM (solid lines) and MEPS (dashed lines) are also shown. 
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Figure 5: The measured scaled momentum distributions (l/N)(n(Kg)/Ax p ) as 
functions of x p in different regions of Q 2 (dots). Other details as in the caption to 
Fig. 4. 
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Figure 6: The measured scaled momentum distributions (l/N)(n(A) / Ax p ) as 
functions of Q 2 in different regions of x p (dots). Other details as in the caption to 
Fig. 4. 
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Figure 7: The measured scaled momentum distributions (l/N)(n(A) / Ax p ) as 
functions of x p in different regions of Q 2 (dots). Other details as in the caption to 
Fig. I 
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;ure 8: The measured scaled momentum distributions (l/N)(n(H)/Ax p ) for 
= Kg (dots), A (triangles) and light charged particles (squares) as functions of 
in different regions of x, p . Other details as in the caption to Fig. 4- 
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